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ABSTRACT 

We present electron collision strengths and their thermally averaged values for the 
forbidden lines of the astronomically abundant doubly- ionized oxygen ion, 0 2+ , in an 
intermediate coupling scheme using the Breit-Pauli relativistic terms as implemented 
in an R-matrix atomic scattering code. We use several atomic targets for the R-matrix 
scattering calculations including one with 72 atomic terms. We also compare with new 
results obtained using the intermediate coupling frame transformation method. We 
find spectroscopically significant differences against a recent Breit-Pauli calculation 
for the excitation of the [O ill] A4363 transition but confirm the results of earlier 
calculations. 

Key words: 0 2+ ion - atomic transition - atomic spectroscopy - collision strength 

- effective collision strength - Upsilon - planetary nebulae - nebular physics - optical 

- infrared - forbidden lines. 


1 INTRODUCTION 


The forbidden lines of 0 2+ are among the most important 
features in the spectra of photoionized plasmas which in- 
clude, inter alia, H il regions and planetary nebulae. The 
exceptional brightness of the strongest [O ill] lines means 
that they can be used to determine the oxygen abundances 
and physical conditions in the Milky Way and other galax- 
ies out to cosmological distances that reach redshifts of more 
than z = 3 I Maiolino et ra/|[2008 ) . 


It has been recently suggested (Nicholls et al 20121 
that the elemental abundance and electron temperature 
anomalies seen in the analysis of the planetary nebula spec- 
tra, where considerable differences have been observed be- 
tween the results obtained from the collisionally-excited 
lines (CEL) and those obtained from the optical recombina- 
tion lines (ORL), might be resolved by using non Maxwell- 
Boltzmann (MB) distributions for the energies of the free 
electrons. The k distribution, which is widely used in the 
analysis of solar data, was proposed as a replacement for the 
MB distribution to resolve this issue. If the electron distribu- 
tions are generally non-Maxwellian in nebulae it would affect 
the analysis of [O ill] lines significantly and reliable collision 
strength data are needed to compute the effective collision 
strengths for collisional excitation and de-excitation. 

The proposal that the electron energy distribution in 
planetary nebulae is not Maxwellian dates back to the 1940s 
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at least where | Hagihara ( 1944j ) proposed that the velocity 
distribution of free electrons in gaseous assemblies, such as 
those found in planetary nebulae, deviates significantly from 
the Maxwellian. Bohm & Al ler| |I94 7 ) argued against Hagi- 
hara and concluded that any deviation from the Maxwellian 
equilibrium distribution is very small. The essence of Bohm 
and Aller’s argument is that for typical planetary nebulae 
conditions of electron temperature of about 10000 K and 
electron number density of about 10 4 cm -3 , the thermal- 
ization process of elastic collisions between an electron and 
other electrons and ions is by far the most frequent event 
and typically occurs once every second, while other pro- 
cesses that shift the system from its thermodynamic equi- 
librium, like inelastic scattering with other ions that leads 
to metastable excitation or recapture, occur at much larger 
time scales estimated to be months or even years. Bohm 
and Aller also indicated the significance of any possible de- 
viation from a Maxwellian distribution on derived elemental 
abundances. 


Although there have been many studies related to col- 
lision strengths of 0 2+ , as we will discuss in the coming 
paragraphs, some of the previous data have limitations. For 
example, some of these data are produced in an LS cou- 
pling scheme while others are based on approaches that do 
not adequately treat resonance phenomena. 

Before the advent of close-coupling codes there were 
several calculations of collision strengths for excitation of 
the O hi forbidden lines that did not incorporate resonance 
effects ( |Czyza k et al|1968| |Seaton|1975| |Bhatia et al|1979| ). 


The first close-coupled collision strengths were obtained 


by Baluja et al (1980) for some of the semi- forbidden inter- 
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combination transitions of O ill using the R-matrix method 
I Berrington et q/|i 97 4 1 |i 9 8 7 [ pTn m m e r et al 1993||Berrington| 


et aZ|| 1995 I. They included all channels with configurations 
Is 2s 2p 2 , Is 2 2s 2p 3 and Is 2 2p 4 in the expansion of the 
wavefunction. They also used three pseudo-orbitals (3s, 3p 
and 3d) and allowed for configuration interaction in the in- 
cluded states with the addition of correlation terms in the 
total wavefunction. 

Ho fe Henry] ( 1983]) also used the close coupling approx- 
imation with configuration interaction in the target wave- 
function to compute the collision strengths of some of O ill 
transitions in LS coupling. They employed a mix of spec- 
troscopic and correlation Hartree-Fock orbitals to describe 
their target. 


Relatively extensive work was done by Aggarwal ( 1983 


1985 ) who computed collision strengths of O ill transitions 
between the fine structure levels using configuration inter- 
action target wavefunctions. He transformed LS coupling 
reactance matrices obtained from R-matrix calculations to 
pair coupling with the program JAJOM ( |Sara ph 1978]. The 
results were obtained with a fine energy mesh up to 5.16 Ry- 
dberg where a complex resonance structure was observed on 
the entire mesh. 


Aggarwal ( 1993 ) used an elaborate configuration inter- 
action target described by Aggarwal & Hibbert ( 1991|) and 


the R-matrix method in LS coupling to compute effective 
collision strengths for some inelastic transitions of O ill be- 
tween 26 LS-coupled states of six configurations over a wide 
range of electron temperature (2500-200000 K). They em- 
ployed the standard and no-exchange R.-matrix codes on a 
fine energy mesh that reveals the resonance structure. This 


work was extended by Aggarwal & Keenan (1999), who they 


computed the collision strengths for the transitions between 
the fine structure levels using the R-matrix method includ- 
ing all partial waves with L ^ 40 to ensure convergence. 


Aggarwal & Keenan 

1 1999 1 transformed the LS reactance 

matrices obtained by 

Aggarwal (1993) into pair coupling us- 

ing JAJOM (Saraph 

1978) where necessary. They only tab- 


ulated fine-structure collision strengths for some transitions, 
pointing out that in pair coupling, if one of the terms in a 
transition has spin zero and hence J = L, e.g. 3 P 4 D, the 
fine-structure collision strengths are proportional to the sta- 
tistical weight of the non-zero spin states, in this example 
the 3 Pj levels. 


Lennon & Burke ( 1994) did extensive work on 0 2+ colli- 
sion strengths for the transitions between the fine structure 
levels, as part of a wider investigation on the carbon iso- 
electronic ions, using the R-matrix method, where the CIV3 


configuration interaction code (Hibbert 19751 was used to 


generate the target wavefunctions. The target included 12 
states belonging to 3 configurations (Is 2 2s 2 2p 2 , Is 2 2s 2p 3 
and Is 2 2p 4 ). They also transformed to pair coupling in the 
same way as 


Aggarwal & Keenan (1999) described above. 


They presented a sample of Maxwellian based effective col- 
lision strengths in the temperature range 10 3 — 10° K. 


Recently, Palay et al |2012 1 made the first calculation of 
collision strengths for the O ill forbidden transitions using a 
relativistic Breit-Pauli (BP) R-matrix method with resolved 
resonance structures. They used 22 configurations (3 spec- 
troscopic and 19 correlation) to describe the target. Like 
most of the previous studies, they have also presented sam- 


ples of the Maxwellian averaged effective collision strengths 
which were also computed at temperatures down to 100 K. 
The most recent R-matrix calculations (|Lennon fc| 


Burke 1994 Aggarwal & Keenan 1999 Palay et al 2012 1 


generally agree to within 10% for the thermally averaged 
collision strengths for the forbidden transitions among the 
five lowest levels. An exception to this generally close agree- 
ment is for the transitions from the lowest three 3 Pj levels 


to the 4 So state. The recent results of Palay et al (2012) 
differ significantly from those of earlier workers. The exci- 
tation mechanism of the 4 So level is important because the 
4 So — t 1 D 2 A4363 line is widely used to infer the electron 
temperature in H II regions and planetary nebulae. If a k 
distribution of electron energies is assumed, the number of 
free electrons capable of exciting the 4 So state would be in- 
creased relative to a MB distribution which would affect the 
derived 0 2+ abundance. 

The aim of the present paper is twofold. Firstly we 
make a Breit-Pauli R-matrix calculation of the 0 2+ colli- 
sion strengths with an independently derived target config- 
uration basis to compare with previous work, especially the 


only other Breit-Pauli results from Palay et al (2012). Sec- 
ondly we attempt to place realistic error estimates on our 
results by examining the effect of several factors on our re- 
sults. We discuss the convergence of our calculation as the 
number of target states is increased. Our largest target in- 
cludes significant contributions to the dipole polarizability 
of the three energetically lowest terms. We also consider the 
effect of Gailitis averaging of the collision strengths close 
to the excitation thresholds, especially for excitation of the 
3 Pi level between the 3 Pi and 3 P 2 thresholds. We addi- 
tionally compare the results of the Breit-Pauli calculation 
with those obtained using the Intermediate Coupling Frame 


Transformation (ICFT) R-matrix method (Griffin, Badnell 
|fe Pin dzola 1998). This method is based on transforming 
the non-physical LS'-coupled reactance matrices, to compute 
collision strengths in intermediate coupling. 

The calculation described in the following sections is 
constructed to provide accurate results for the excitation of 
the optical and infrared forbidden transitions among the five 
lowest levels of O a+ at temperatures typical of photoionized 
plasmas. We compute collision strengths up to « 1.3 Ry- 
dberg free electron energy relative to the ground level and 
Mawell-Boltzmann averaged collision strengths from 100 K 
to 25000 K. 

The main too ls used in this inv estigation are the Au- 
tostructure codJ^ (Eissner et a(|l974 Nussbaumer fe Storey| 


19781 Badnell 2011) to define and elaborate the atomic 
target and the UCL-Belfast-Strathclyde R-matrix codtj^] 
I Be rrington et a/||1995 I to do the actual scattering calcula- 
tions. We compare our results with earlier calculations and 
also assess the reliability of our results. 


1 See Badnell: Autostructure write-up on WWW. URL: 

amdpp.phys. strath. ac.uk/autos/ver/WRITEUP. 

2 See Badnell: R-matrix write-up on WWW. URL: 

amdpp.phys.strath.ac.uk/UK_RmaX/codes/serial/WRITEUP. 
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Table 1. The configuration basis used to define the scattering 
target. The Is 2 core is to be understood in all configurations. 
The bar signifies a correlation orbital. 


2s 2 2p 2 , 2s 2p 3 , 2p 4 

2s 2 2p 31; l = 0, 1.2 

2s 2p 2 31; 1 = 0,1,2 

2p 3 31; 1 = 0, 1,2 

2s 2 31 3 1', l,l' =0,1,2 

2s 2p 31 3 1', 1, l' =0,1,2 

2p 2 31 3 1', 1,/' =0,1,2 

2s 2 2p 4f, 2s 2p 2 4f, 2p 3 4f 

2s 2 3d 4f, 2s 2p 3d 4f, 2p 2 3d 4f 

2s 2 If 2 , 2s 2p 4f 2 , 2p 2 4f 2 


Table 2. Orbital scaling parameters, A n ;, for Autostructure in- 
put. The rows stand for the principal quantum number n, while 
the columns stand for the orbital angular momentum quantum 
number 1. 



s 

P 

d 

f 

1 

1.44889 




2 

1.22418 

1.18282 



3 

-0.80508 

-0.61905 

-1.04731 


4 




-1.87410 


2 COMPUTATION 

In the following we outline the computational methods used 
in this work. 


2.1 The 0 2+ target 

We used the Autostructure code (Badnell 2011 1 to generate 
the target radial functions required as an input to the first 
stage of the R-matrix code. The radial data were generated 
using thirty-nine configurations containing seven orbitals; 
three physical (Is, 2s and 2p) and four correlation orbitals 
(3s, 3p, 3d, 4f). These configurations are given in Table [l] 
An iterative optimization variational protocol was used to 
obtain the orbital scaling parameters, A „i, which are given in 
Table[2] The correlation orbitals are calculated in a Coulomb 
potential with central charge 8|A n !|. 

In the scattering calculations, targets with differing 
numbers of target states were used, with the largest hav- 
ing 72 terms which are listed in Table [3] Calculations were 
also made with the first 10 and 20 terms from this list as 
discussed in more detail below. 

In Table [4] we show the statistically weighted oscilla- 
tor strengths, gf, in the length and velocity formulations 
for all the strong allowed transitions between the 2s 2 2p 2 
and 2s 2p 3 configurations. The agreement is excellent. Good 
agreement between the length and velocity results is a nec- 
essary but not sufficient condition for ensuring the quality 
of the target wave functions. These transitions also make 
the largest contributions to the dipole polarizabilities of the 
three lowest terms. 

The 72 terms listed in Table [3] were chosen to include 
all those correlation configurations that contribute signifi- 
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cantly to the dipole polarizability of the three lowest terms. 
The main contributions come from 2s 2 2p 3d configuration. 
The contribution of states outside the n — 2 complex to the 
dipole polarizabilities of the 3 P, : D and 3 S terms is 37%, 
37% and 60% respectively. In Table [5] we list the energies 
of the 18 levels of the n = 2 complex configurations. We 
show theoretical energies which include one- and two-body 
fine-structure interactions (E t h 2 ) and those which only in- 
clude the spin-orbit interaction (Sthi), the latter being the 
only fine-structure interactions included in the version of the 
R-matrix code that we use (see footnote 2). We return to 
the effect of omitting two-body fine structure interactions in 
section |3] 


2.2 The Scattering Calculations 


We made several calculations with increasing numbers of 
target states, both with Breit-Pauli and the Intermediate 
Coupling Frame Transformation R-matrix methods. The 
target radial functions were supplied as a radial grid for- 
mat rather than Slater type orbital format where the radial 
file was generated by Autostructure. The inner region ra- 
dius (RA) in the R-matrix formulation was 9.315 au. Twelve 
continuum basis functions were used to represent the wave- 
functions in the inner region. This choice was based on con- 
vergence tests and with experience from previous work on 


the C + +e system (Sochi 


2012 


Sochi & Storey 2013). The 


maximum value of 2J for the ( N + 1) -electron problem was 
chosen to be 19 although 15 was found to be sufficient for 
convergence of the collision strengths for the forbidden tran- 
sitions of interest here. 

As indicated previously, we made three sets of calcula- 
tions using the configuration basis described in section [2 T] 
with 10- , 20- and 72-terms using both the BP and ICFT 
approaches. These three targets comprise 18, 34 and 146 
fine-structure levels respectively. The main purpose of using 
several targets is to have an estimate of the error in the final 
results from observing the convergence of the results with 
different numbers of target terms. For all three targets, the 
( N + l)-electron wavefunction contains all possible configu- 
rations formed from the 39 configurations of the A-electron 
target combined with any of the orbitals, spectroscopic and 
correlation. There are 102 such ( N + l)-electron configura- 
tions. 

Experimental energies obtained from the National In- 
stitute of Standards and Technology (NIST([^] were used in 
place of theoretical ones to ensure correct positioning of 
thresholds for convergence of resonance series. In some cases 
this required re-ordering the target states. 

Collision strengths were calculated for electron energies 
up to 1.28 Rydberg relative to the 2s 2 2p 2 3 P° ground level, 
hence 0.89 Rydberg relative to the highest state of interest, 
2s 2 2p 2 1 S®. This energy corresponds to ~ 7 kT when the 
electron temperature T = 20000 K, the approximate upper 
limit for temperatures in photoionized nebulae. Over this 
energy range, collision strengths were calculated at 20000 
equally spaced energies, except between the 2s 2 2p 2 3 P® and 
3 P® levels where calculations were performed on a mesh 100 
times finer. Calculations were also made with and without 


3 See NIST website: www.nist.gov. 
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Table 3. Target terms and energies, E , calculated by Autostructure using the configuration basis listed in Table [T] The Is 2 core is 
suppressed from all configurations. All these terms are included in the 72-term target, while for the smaller targets (10- and 20-term) 
only the first 10 and 20 terms respectively are included. 


Index 

Configuration 

Term 

E (cm 1 

Index 

Configuration 

Term 

E (cm 3 ) 


1 

2s 2 2p 2 

3pe 

0.0 

37 

2s 2p 2 3p 

1 s° 

483995 


2 

2s 2 2p 2 

1 D° 

21257 

38 

2s 2p 2 3p 

3po 

485363 


3 

2s 2 2p 2 

i S e 

45630 

39 

2s 2p 2 3s 

3 S e 

486363 


4 

2s 2p 3 

5 S° 

58948 

40 

2s 2p 2 3p 

lpo 

486594 


5 

2s 2p 3 

3 d° 

121133 

41 

2s 2p 2 3p 

3 d° 

491627 


6 

2s 2p 3 

3po 

144640 

42 

2s 2p 2 3s 

3pe 

496357 


7 

2s 2p 3 

1 D° 

190314 

43 

2s 2p 2 3p 

3po 

498566 


8 

2s 2p 3 

3 S° 

199693 

44 

2s 2p 2 3p 

3 s° 

503468 


9 

2s 2p 3 

1 po 

214885 

45 

2s 2p 2 3s 

1 pc 

506149 


10 

2 P 4 

3pe 

287613 

46 

2s 2p 2 3s 

ige 

508564 

11 

2s 2 2p 3p 

lpe 

301182 

47 

2s 2p 2 3p 

1 D° 

509665 

12 

2 P 4 

1 D e 

302782 

48 

2s 2p 2 3p 

lpo 

526999 

13 

2s 2 2p 3p 

3 D h 

306265 

49 

2s 2 2p 3d 

3po 

539361 

14 

2s 2 2p 3s 

3po 

309248 

50 

2p 3 3p 

5pe 

542552 

15 

2s 2 2p 3p 

3ge 

310788 

51 

2s 2 2p 3d 

1 D° 

550865 

16 

2s 2 2p 3p 

3pe 

315730 

52 

2p 3 3s 

5 S° 

551832 

17 

2s 2 2p 3s 

1 po 

323156 

53 

2p 3 3p 

3pe 

553504 

18 

2s 2 2p 3p 

1 D B 

328464 

54 

2p 3 3p 

3 D e 

567272 

19 

2s 2 2p 3p 

i S e 

345567 

55 

2p 3 3p 

lpe 

567946 

20 

2p 4 

4 S e 

351203 

56 

2p 3 3p 

3pe 

568179 

21 

2s 2p 2 3p 

3 S° 

372370 

57 

2p 3 3p 

lpe 

570684 

22 

2s 2p 2 3p 

5 d° 

376169 

58 

2s 2 2p 3d 

3po 

572129 

23 

2s 2p 2 3s 

5pe 

379161 

59 

2s 2 2p 3d 

3 d° 

578854 

24 

2s 2p 2 3p 

5po 

379976 

60 

2p 3 3s 

3 d° 

584855 

25 

2s 2p 2 3p 

3 d° 

392038 

61 

2p 3 3s 

3g° 

589631 

26 

2s 2p 2 3p 

5 s° 

395094 

62 

2p 3 3p 

3pe 

600736 

27 

2s 2p 2 3p 

3po 

400796 

63 

2s 2p 2 3d 

5pe 

600824 

28 

2s 2p 2 3s 

3pe 

418063 

64 

2p 3 3s 

1 D° 

602294 

29 

2s 2p 2 3p 

3po 

437947 

65 

2p 3 3p 

1 D e 

606079 

30 

2s 2p 2 3p 

1 D° 

439915 

66 

2p 3 3p 

Sge 

606299 

31 

2s 2p 2 3s 

3 D e 

442045 

67 

2p 3 3p 

3 D e 

607292 

32 

2s 2p 2 3p 

lpo 

443431 

68 

2s 2p 2 3d 

5 D e 

610073 

33 

2s 2p 2 3p 

3 D° 

447068 

69 

2p 3 3p 

lpe 

612417 

34 

2s 2p 2 3p 

1 po 

449383 

70 

2s 2 2p 3d 

lpo 

617417 

35 

2s 2p 2 3p 

3po 

456397 

71 

2p 3 3p 

3pe 

618955 

36 

2s 2p 2 3s 

1 D e 

468076 

72 

2s 2 2p 3d 

lpo 

620437 


I ^ ^ xj | f 1/2 

Table 4. The weighted LS oscillator strengths in length, (gf) L , ( 6j , T f ) = — eVW ) J°° n y ( e< ) ( — ) / © , T f ) ds (1) 

and velocity, ( gf ) v , forms for the transitions within the n = 2 2 0 

complex. 


(s/)i (ff/)v 


3pe 

!D c 

igo 


3 D ° 

3po 

3 s° 

iD" 

lpo 

lpo 


1.01 

1.28 

1.70 

1.57 

1.18 

0.27 


0.99 

1.31 

1.61 

1.61 

1.22 

0.27 


where Tf is the effective temperature, k is the Boltz- 
mann constant, ti and tj are the free electron energy rela- 
tive to the states i and j respectively, A Eij (= tj — ti ) is the 
energy difference between the two states, 0,:j is the collision 
strength of the transition between the i and j states, and 
f{ti, Tf) is the energy- and temperature-dependent electron 
distribution. In what follows we will only consider Maxwell- 
Boltzmann distributions of electron energy, given by 


Gailitis averaging of the collision strengths in the region be- 
neath each threshold where the effective quantum number 
v > 10. 

We calculate the thermodynamically-averaged collision 
strengths for electron excitation, T, from a lower state i to 
an upper state j from 


/mb(£,T) 


2 

(. fcT ) 3/ 2 



kT 

5 


(2) 


although, as discussed above, other distributions such as the 


k distribution have been proposed and discussed ( Vasyliunas 
1968 Nicholls et a/|2012 Storey fc Sochi|2013 l. 
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(a) 1-2 


(b) 1-3 


(c) 1-4 









Figure 2. Thermally averaged collision strengths from the 10-term target (solid line) and 20-term target (dotted line), shown as the 
percentage difference from the 72-term target, plotted against log T [K]. The labels of the sub-figures refer to the level indices in Table 

m 


3 RESULTS AND DISCUSSION 
3.1 Results 

A sample of our 10-, 20- and 72-term BP collision strengths 
is shown in Figure [l] The agreement between the three cal- 
culations is excellent, with the most obvious difference being 
that some resonances move to lower energies as the target 
size is increased, as might be expected. Figure [2] shows the 
results for the thermally averaged collision strength, T, as a 
function of temperature for the 10- and 20-term calculations 
relative to the 72-term calculation as a percentage difference. 
The differences are less than 9% at any temperature for the 
10-term calculation and less than 5% for the 20-term case. 


The energy region between the 2s 2 2p 2 d P° and 3 P“ 
states contains a Rydberg series of resonances converging 
on the 3 P® level with an effective quantum number at the 
3 P® threshold of 47.7. The energy difference between the 3 P° 
and 3 P® levels of 193 cm -1 corresponds to a temperature of 
278 K, so this energy region is significant for computing T 
at temperatures down to 100 K. We calculate the collision 
strengths with an energy interval of 6.4 x 10~' Rydberg in 
this interval and compare with the result of using Gailitis 
averaging in this region. The difference is less than 1% at 
any temperature and we conclude that Gailitis averaging is 
adequate to obtain accurate values of T down to 100 K. 

The results of the ICFT calculations showed unexpect- 


© 0000 RAS, MNRAS 000, 000-000 


6 P.J. Storey, Taha Sochi & N.R. Badnell 


Table 5. The 18 lowest energy levels of 0 2+ and their experimen- 
tal (E e x) and theoretical (Ethi and E t ^o) energies in wavenum- 
bers (cm -1 ). The experimental energies are obtained from the 
NIST database while the theoretical energies were obtained from 
Autostructure with the configuration basis listed in Table [T] The 
energies E t hi were obtained with only spin-orbit terms in the tar- 
get Hamiltonian while -Eth2 also include two-body fine-structure 
interactions within the n = 2 complex. 


Index 

Level 

L'ex 

Ethl 

®th2 

1 

2s 2 2p 2 3 P® 

0.00 

0 

0 

2 

2s 2 2p 2 3 P® 

113.18 

115 

113 

3 

2s 2 2p 2 3 P® 

306.17 

339 

308 

4 

2s 2 2p 2 

20273.27 

21489 

21471 

5 

2s 2 2p 2 iS® 

43185.74 

45900 

45882 

6 

2s 2p 3 5 S° 

60324.79 

59600 

59582 

7 

2s 2p 3 3 D° 

120058.2 

121800 

121775 

8 

2s 2p 3 3 D° 

120053.4 

121804 

121799 

9 

2s 2p 3 3 D° 

120025.2 

121812 

121805 

10 

2s 2p 3 3 P° 

142393.5 

145316 

145307 

11 

2s 2p 3 3 P° 

142381.8 

145321 

145307 

12 

2s 2p 3 3 P° 

142381.0 

145331 

145319 

13 

2s 2p 3 4 D° 

187054.0 

191025 

191006 

14 

2s 2p 3 3 S° 

197087.7 

200423 

200405 

15 

2s 2p 3 ip° 

210461.8 

215609 

215591 

16 

2p 4 3 P® 

283759.7 

288653 

288629 

17 

2p 4 3 P® 

283977.4 

288863 

288854 

18 

2p 4 3 P® 

284071.9 

288965 

288951 


edly large differences from the BP results in some energy 
domains. This is illustrated in Figure [3] where we compare 
the thermally averaged collision strengths for the 72-term 
ICFT calculations with the 72-term BP results for the 3 P® 
3 P® transition. Due to the difference in scaling with effec- 
tive charge (z e ff) of term energy separations (oc z e g) and 
resonance energies (oc z^g ) , resonance effective quantum 
numbers can become small for lowly ionized systems. Such 
deeply-closed channels can be problematic for the multi- 
channel quantum defect theory (MQDT) used by the ICFT 
method due to computational finite numerical precision of 
highly divergent wavefunctions. Gorczyca & Badnell] (2000 1 
found that classically forbidden channels (e.g. n < l) could 
be handled expediently by simply omitting them from the 
MQDT representation. For low-energy scattering in O a+ we 
encountered a similar problem in a new guise for n < 2. 
The closed-channel partition of the MQDT representation 
should give no contribution since all bound orbitals (spec- 
troscopic and pseudo) are projected out of the continuum 
basis. All such closed channel contributions (e.g. correla- 
tion resonances) arise instead in the open-open part of the 
scattering matrix. For l > 1 the original Gorczyca fe Bad-| 
|nell | ( |2000[ > expediency already omits such closed channels 
(n < l). For l = 0, 1 we found it necessary to explicitly omit 
such closed channels from the closed partition as well. We 
show the effect of this modification as the dashed line in Fig- 
ure [3] The agreement with the full Breit-Pauli calculation is 
now excellent. 

Considering the convergence as the number of target 
states is increased and the good agreement between the 
ICFT and Breit-Pauli results, we adopt the results of the 72- 
term Breit-Pauli calculation as our final results and, based 





(c) 1-4 

Figure 1. Collision strength (vertical axis) versus final electron 
energy in Rydberg (horizontal axis) for Breit-Pauli calculations 
of the 1-2, 1-3 and 1-4 transitions with 10-term (dotted), 20-term 
(dashed), and 72-term (solid). Refer to Table[5]for level indexing. 


on the convergence behavior and the effect of Gailitis aver- 
aging, estimate an uncertainty of no more than 5% in the 
final thermally averaged collision strengths. In Table [6] we 
tabulate thermally averaged collision strengths T, for the 72- 
term target in the temperature range logioT = 2.0(0. 1)4.4. 


3.2 Comparison to Previous Work 


We compare our effective collision strength results with 
those from previous calculations of similar quality, that is 
those which used close-coupling techniques and computed 
collision strengths at sufficient energies to delineate reso- 
nances. 

In Table [7] we compare our final 72-term results with 
the LS results of |Lennon fe Burke (19941. That calculation 
was based on the 12-state target including n = 3 correlation 


orbitals described by Burke, Lennon & Seaton (1989). They 


agree within 10% for all transitions and all temperatures. 
The effective collision strengths, T( 3 P - 2 D) and X( 3 P - 
1 S), for excitation of the optical forbidden lines do not differ 
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Table 6. Thermally averaged collision strengths from the 72-term Breit-Pauli calculation as a function of temperature. The transitions 
are indexed as in Table [5] 


log T [K] 

1-2 

1-3 

1-4 

1-5 

2-3 

2-4 

2-5 

3-4 

3-5 

4-5 

2.0 

0.63500 

0.22586 

0.23177 

0.02989 

1.11212 

0.69748 

0.09003 

1.17021 

0.15117 

0.38267 

2.1 

0.62649 

0.22561 

0.23211 

0.02989 

1.10993 

0.69849 

0.09003 

1.17188 

0.15116 

0.38291 

2.2 

0.61540 

0.22529 

0.23259 

0.02989 

1.10711 

0.69994 

0.09001 

1.17427 

0.15112 

0.38316 

2.3 

0.60216 

0.22490 

0.23335 

0.02988 

1.10362 

0.70220 

0.08998 

1.17802 

0.15106 

0.38346 

2.4 

0.58741 

0.22447 

0.23449 

0.02987 

1.09952 

0.70562 

0.08994 

1.18370 

0.15099 

0.38382 

2.5 

0.57193 

0.22411 

0.23600 

0.02985 

1.09530 

0.71013 

0.08989 

1.19123 

0.15089 

0.38428 

2.6 

0.55668 

0.22408 

0.23761 

0.02983 

1.09235 

0.71497 

0.08982 

1.19932 

0.15075 

0.38486 

2.7 

0.54292 

0.22475 

0.23888 

0.02980 

1.09280 

0.71876 

0.08972 

1.20571 

0.15058 

0.38561 

2.8 

0.53194 

0.22633 

0.23934 

0.02976 

1.09839 

0.72014 

0.08960 

1.20807 

0.15036 

0.38660 

2.9 

0.52440 

0.22869 

0.23872 

0.02970 

1.10890 

0.71829 

0.08942 

1.20503 

0.15007 

0.38792 

3.0 

0.51992 

0.23126 

0.23701 

0.02963 

1.12183 

0.71317 

0.08919 

1.19654 

0.14969 

0.38974 

3.1 

0.51730 

0.23343 

0.23441 

0.02952 

1.13386 

0.70539 

0.08887 

1.18365 

0.14917 

0.39242 

3.2 

0.51544 

0.23489 

0.23123 

0.02937 

1.14303 

0.69591 

0.08842 

1.16800 

0.14846 

0.39684 

3.3 

0.51396 

0.23583 

0.22786 

0.02918 

1.14971 

0.68588 

0.08787 

1.15154 

0.14756 

0.40497 

3.4 

0.51320 

0.23671 

0.22485 

0.02900 

1.15583 

0.67700 

0.08732 

1.13712 

0.14670 

0.41997 

3.5 

0.51369 

0.23794 

0.22306 

0.02891 

1.16333 

0.67179 

0.08709 

1.12887 

0.14636 

0.44484 

3.6 

0.51585 

0.23978 

0.22342 

0.02906 

1.17355 

0.67301 

0.08756 

1.13139 

0.14719 

0.47990 

3.7 

0.51993 

0.24246 

0.22653 

0.02954 

1.18745 

0.68249 

0.08901 

1.14765 

0.14967 

0.52127 

3.8 

0.52604 

0.24618 

0.23236 

0.03036 

1.20566 

0.70010 

0.09148 

1.17740 

0.15386 

0.56212 

3.9 

0.53378 

0.25101 

0.24025 

0.03144 

1.22793 

0.72382 

0.09475 

1.21720 

0.15940 

0.59558 

4.0 

0.54210 

0.25678 

0.24917 

0.03268 

1.25256 

0.75062 

0.09849 

1.26198 

0.16570 

0.61740 

4.1 

0.54956 

0.26315 

0.25817 

0.03395 

1.27703 

0.77758 

0.10233 

1.30687 

0.17219 

0.62672 

4.2 

0.55511 

0.26978 

0.26652 

0.03517 

1.29937 

0.80262 

0.10601 

1.34841 

0.17837 

0.62540 

4.3 

0.55866 

0.27654 

0.27384 

0.03625 

1.31924 

0.82453 

0.10927 

1.38465 

0.18385 

0.61639 

4.4 

0.56091 

0.28349 

0.27983 

0.03710 

1.33776 

0.84246 

0.11183 

1.41423 

0.18816 

0.60225 



T[K] 


Figure 3. Effective collision strength (Upsilon) versus tempera- 
ture from the 72-term target: Breit-Pauli (solid line), ICFT before 
modification (dotted line) and ICFT after modification (dashed 
line), for the transition between levels 2 and 3 of Table[5] On this 
scale the solid and dashed lines are almost indistinguishable. 


by more than 6% at any temperature. The agreement is gen- 
erally even better with our 20-term calculation which might 
be expected since that calculation includes the 12 terms of 
the n = 2 complex which is the target of |Lennon fe Burke| 
( |1994[ ). However, their target does not include the states 
constructed from correlation orbitals that make a large con- 
tribution to the polarizability of the important states, as 
discussed in section [2~TI 

The most recent R-matrix calculations where fine- 
structure collision strengths are presented are those of |Ag-| 


garwal & Keenan (1999) and Palay et al (20121. The for- 


mer calculation is based on an elaborate 26-term target de- 


scribed by Aggarwal & Hibbert ( 1991 1 constructed from Is, 
2s and 2p spectroscopic and 3s, 3p, 3d, 4s, 4p and 4d cor- 
relation orbitals. The resulting LS coupled reactance matri- 


ces were recoupled algebraically using the JAJOM (Sarapli 


1978 1 program where necessary. This approach neglects the 


fine-structure interactions between target states and in this 
approximation some fine-structure collision strengths can be 
derived directly from LS-coupled collision strengths using 
only statistical weight factors as described by bot h |Aggar- 
wal & Keenan ( 1999 1 and Lennon & Burke ( 1994 1 . Palay et 


al (2012) have made a 19-level Breit-Pauli R.-matrix calcu- 


lation where the target is expanded over a configuration set 
involving Is, 2s, 2p and 3p spectroscopic orbitals and 3p, 3d, 
4s and 4p correlation orbitals. P alay et~ai\ (2012) use an ex- 
tended version of the Breit-Pauli R-matrix code which they 
attribute to Eissner & Chen (in preparation) which includes 
two-body fine-structure interactions which enables them to 
calculate the fine structure splitting of the ground 3 Pj lev- 
els with an error of order 3% ( Palay et al |2012| ). [Pal ay et al\ 
(20121 were also the first to extend the tabulation of ther- 


mally averaged collision strengths down to very low electron 
temperatures (100 K). 

In Figure [4] we compare graphically our fine-structure 
results with those of Lennon & Burke ( 1994 1 , Aggarwal & 
Keenan ( 1999| ) and |Palay et al l 2012 1 . In Table|8]we compare 
the same results numerically and also include the results of 


the earlier R-matrix calculation by Aggarwal (1983). Fig- 


ure [4] shows the percentage difference in the thermally av- 
eraged collision strengths from these three calculations rel- 
ative to our results, for all ten transitions among the en- 
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ergetically lowest five levels. Where necessary, we derived 
fine-structure collision strengths from the results of|L ennon 


& Burke (19941 and Aggarwal & Keenan (1999) using sta- 


tistical weight factors as outlined above. With the exception 
of the - 1 S|( t ransition, our results agree with those of 


Lennon & Burke 1994 ) and Aggarwal & Keenan (1999) 


to within 10% for all temperatures between 1000 K and 
25000 K where comparison can be made and to within 5% 
for the majority of temperatures. For these two calculations 
the differences are relatively insensitive to temperature, in- 
dicating that their collision strengths have a similar energy 
dependence to ours. We find generally larger disagreements 
with the results of Palay et al (2012), reaching 10-15% at 
the extremes of tabulated temperature for many transitions 
and being even larger for the transitions from the ground 
3 Pj levels to the x So state (transitions 1-5, 2-5 and 3-5). 
Here the differences reach 100% at 100 K and are over 20% 
at 10000 K. The differences also show a distinctive tempera- 
ture dependence. With the exception of the 3 Pj - 3 So tran- 
sitions the Palay et al (2012) results are generally smaller 
than ours at the lowest temperatures and larger at the high- 
est temperatures. This suggests that their collision strengths 
generally have a different energy dependence in the energy 
range relevant for nebular temperatures, about 1/4 Rydberg 
above threshold. 


3.3 Discussion 


In photoionized plasmas the O ill forbidden lines are com- 
monly used to determine the electron temperature of the 
emitting material, and hence to determine the number of 
0 2+ emitters relative to H by comparison with a strong H 
recombination line. The temperature determination rests on 
the ratio of the intensity of the A4363 line to either or both 
of the A4959 and A5007 lines. The A4363 line is relatively 
weak and cannot be seen if the temperature is much below 
5000 K. Once the temperature is known, the much stronger 
AA4959, 5007 lines can be used to deduce the 0 2+ number 
density. In nebular plasmas all these lines are excited colli- 
sionally from the 3 Pj ground levels. The excitation mecha- 
nism for A4363 is therefore central to determining the elec- 
tron temperature and abundances. In Figure [5] we show how 
the derived electron temperature from our work differs from 
that obtained from Lennon fe Burke| ( |1994|) and fr om the 
data of Aggarwal & Keenan (19991 and Palay et al (2012 1 . 


In all the temperature determinations the radiative tran- 
sition probabilities were taken from |Nu ssbaume r fe Storeyl 
( 1981 1 and Storey & Zeippen ( 20001. Very similar tempera- 


tures are obtained with the collision strength data of|Lennon 


& Burke (19941 and Aggarwal & Keenan (19991. Palay et 
~al f ]2012|) state that there are no significant differences in 


line ratios arising from their calculation when comparing to 
Aggarwal & Keenan (19991 but Figure [5] shows that this is 


not the case. The difference in derived temperature is 213 K 
at 5000 K, 421 K at 10000 K and 504 K at 15000 K. 

In summary, our new Breit-Pauli R-matrix calculation 
generally shows much better agreement for thermally av- 
eraged collision strength with the earlier non- Breit-Pauli R- 


matrix calculations of Lennon & Burke ( 1994 1 and Aggarwal 


& Keenan (1999) than the more recent Breit-Pauli work of 


Palay et al (2012). The results of the important forbidden 


line diagnostic line ratios show the same pattern. The rea- 



T [K] 


Figure 5. The difference in derived electron temperature from 
the A4363/(A4959+A5007) line intensity ratio using the data of 
|Lennon &; B urke ( 199 4]l (solid l ine) , | Aggarwal fe Keenan ( |1999[ > 
(dashed line) and |Palay et aJ\ ( |2012[ > (dotted line) against the 
temperature derived from the present results. 


sons why the |Palay et <zZ| ([2012 ) results differ is not clear. One 
question that arises is whether the two-body fine-structure 
terms that are included in the Breit-Pauli R-matrix formu- 
lation of|Pal ay et ( [2012[ ) and not in our calculation might 
be the cause. We do not believe that this is the case for the 
following reason. In Figure [6] we show two sets of results for 
the thermally averaged collision strength for the 3 P® - 3 P® 
transition from 72-term ICFT calculations. The solid line 
includes the effects of the spin-orbit interaction in the tar- 
get, introduced via the so-called Term-Coupling Coefficients 
(TCCs), while the dashed line shows the results obtained in 
pair-coupling, i.e. without TCCs. Except at the lowest tem- 
peratures ( T < 300 K) they differ by no more than 1%. The 
larger difference at the lowest temperatures simply reflects 
the fact that the pair-coupling calculation does not separate 
the 3 P j levels in energy and therefore the threshold energies 
of these levels are not correct. The results for the other tran- 
sitions show similar behavior. We emphasize, however, that 
the ICFT calculation which does incorporate target spin- 
orbit effects agrees with the full Breit-Pauli calculation to 
within 1% at all temperatures. The good agreement that we 
find shows that the spin-orbit interaction has a very small 
effect on the results. In 0 2+ two-body fine-structure inter- 
actions are substantially smaller than the spin-orbit inter- 
action and should therefore have a negligible effect on the 
results. This point is emphasized in Figure [7] where we show 
the percentage difference between the results of jPalay et al\ 
12012) and ours for the three 3 Pj - 3 So transitions. Except 
at very low temperatures, they do not show any significant 
dependence on J which might be expected if fine-structure 
effects were important and indicate rather that the term- 
term 3 P - 3 S collision strengths differ significantly between 
the two calculations. 


4 CONCLUSIONS 

In the present paper, the collision strengths for the transi- 
tions between the lowest five levels of the astronomically- 
important 0 2+ + e~ atomic system up to about 1.3 Ryd- 
berg of electron excitation energy are computed in the close 
coupling approximation using the UCL-Belfast-Strathclyde 
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Figure 4. Percentage differences of thermally averaged collision strengths from our 72-term Breit-Pauli calculation, (vertical axis) versus 
temperature in 10000 K (horizontal axis). Results from [Lennon &:~ Burkc ( 1994|l (solid line), | Aggarwal K' Keenan (1999) (dashed line) 
and |Palay et al\ (|2012[) (dotted line) . The labels of the sub-figures refer to the level indices in Table |5| 


R-matrix atomic code. Different coupling schemes with dif- 
ferent atomic definitions and parameters are used to describe 
the scattering target and scattering process. 


Our results were extensively compared to previous 
work. We found a good agreement in most cases which in- 
creases our confidence in our results. However, we found sig- 
nificant differences with Palay et al (2012) who also used a 
Breit-Pauli coupling scheme and hence a better agreement 
was expected. The good agreement between our R-matrix 
Breit-Pauli calculation and earlier R-matrix work in which 
the fine-structure was treated more approximately strongly 
supports our results. We showed that the relatively large dif- 


ferences found for the excitation of the A4363 line between 
the work of Palay et al ( 2012 1 on the one hand, and all previ- 
ous calculations, on the other, leads to significant differences 
in derived temperatures from the main [O ill] line ratios. 


With regard to the use of the ICFT method, for lowly 
ionized systems some resonances can have very low princi- 
pal quantum number, and channels are deeply closed, which 
can cause problems for multichannel quantum defect theory. 
This difficulty can be overcome by explicitly omitting chan- 
nels with very low effective quantum number and in any case 
evaporates as the effective charge number increases. 
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Figure 6. Effective collision strength (Upsilon) versus temper- 
ature from 72-term target: ICFT with spin-orbit interactions in 
the target (solid line) and ICFT in pair coupling (dashed line) for 
the transition between levels 2 and 3 of Table [5j 



Figure 7. Percent difference in T versus logarithm of temper- 
ature between the results of P alay et al\ ( |2012| l and our 72- 
term Breit-Pauli calculation for the transitions 3 P®- 1 S® (trian- 
gle), 3 P®- 1 S® (circle) and 3 P“- 1 S® (square). 
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Table 7. Comparison of thermally averaged collision strengths, T, from [Lennon &; Burke| ( |1994| ) (LB) and the current work within the 
lowest five levels as a function of temperature. The first row of each temperature is from Table 3 of LB and the second row is from the 
current work using the Breit-Pauli 72-term target. 


log T [K] 

3pe 3pe 
0 1 

3pe 3pe 
0 2 

3pe 3pe 
1 2 

3pe_lj~)e 

3 P e - 1 S® 

lne Ice 

‘V s o 

3.0 

0.4975 

0.2455 

1.1730 

2.2233 

0.2754 

0.4241 


0.5199 

0.2313 

1.1218 

2.1331 

0.2667 

0.3897 

3.2 

0.5066 

0.2493 

1.1930 

2.1888 

0.2738 

0.4268 


0.5154 

0.2349 

1.1430 

2.0811 

0.2643 

0.3968 

3.4 

0.5115 

0.2509 

1.2030 

2.1416 

0.2713 

0.4357 


0.5132 

0.2367 

1.1558 

2.0237 

0.2610 

0.4200 

3.6 

0.5180 

0.2541 

1.2180 

2.1117 

0.2693 

0.4652 


0.5158 

0.2398 

1.1736 

2.0107 

0.2616 

0.4799 

3.8 

0.5296 

0.2609 

1.2480 

2.1578 

0.2747 

0.5232 


0.5260 

0.2462 

1.2057 

2.0913 

0.2732 

0.5621 

4.0 

0.5454 

0.2713 

1.2910 

2.2892 

0.2925 

0.5815 


0.5421 

0.2568 

1.2526 

2.2425 

0.2941 

0.6174 

4.2 

0.5590 

0.2832 

1.3350 

2.4497 

0.3174 

0.6100 


0.5551 

0.2698 

1.2994 

2.3987 

0.3165 

0.6254 

4.4 

0.5678 

0.2955 

1.3730 

2.5851 

0.3405 

0.6090 


0.5609 

0.2835 

1.3378 

2.5184 

0.3339 

0.6022 
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Table 8. Comparison of thermally averaged collision strengths, T, between |AggarwaI| ( |1983| > (A), |Lennon fc Bur ke (1994) (LB), [Aggarwal | 
|fc Keenan| ( |1999[ > (AK) , |Palay et d\ ( |2012~| (P), and the current work (SSB) using the 72-term target as a function of temperature [K] . 
The transitions are indexed as in Table |5 Note that the values attributed to |Lennon &; B urke (1994) for T = 2500K and 25000K are 
those tabulated for log T = 3.4, 4.4 in that work. 


Index 






Temperature [K] 


100 

500 

1000 

2500 

5000 

7500 10000 12500 15000 17500 20000 25000 30000 


1-2 

A 




0.5041 

0.5172 

0.5310 

0.5417 

0.5490 

0.5537 

0.5567 

0.5586 

0.5612 

0.5633 


LB 



0.4975 

0.5115 



0.5454 





0.5678 



AK 




0.5011 

0.5084 

0.5159 

0.5222 

0.5266 

0.5294 

0.5311 

0.5324 

0.5348 

0.5380 


P 

0.5814 

0.5005 

0.4866 


0.5240 


0.5648 




0.6007 


0.6116 


SSB 

0.6350 

0.5430 

0.5199 

0.5132 

0.5199 

0.5317 

0.5421 

0.5494 

0.5540 

0.5569 

0.5587 

0.5609 

0.5623 

1-3 

A 




0.2499 

0.2566 

0.2646 

0.2717 

0.2776 

0.2824 

0.2865 

0.2901 

0.2962 

0.3013 


LB 



0.2455 

0.2509 



0.2713 





0.2955 



AK 




0.2406 

0.2449 

0.2512 

0.2573 

0.2626 

0.2669 

0.2707 

0.2739 

0.2798 

0.2855 


P 

0.2142 

0.2153 

0.2234 


0.2469 


0.2766 




0.3106 


0.3264 


SSB 

0.2259 

0.2247 

0.2313 

0.2367 

0.2424 

0.2497 

0.2568 

0.2629 

0.2682 

0.2727 

0.2766 

0.2833 

0.2890 

1-4 

A 




0.2283 

0.2262 

0.2337 

0.2426 

0.2506 

0.2627 

0.2627 

0.2672 

0.2740 

0.2790 


LB 



0.2470 

0.2380 



0.2544 





0.2872 



AK 




0.2260 

0.2265 

0.2343 

0.2434 

0.2515 

0.2582 

0.2637 

0.2683 

0.2751 

0.2799 


P 

0.1959 

0.2088 

0.2154 


0.2347 


0.2693 




0.3094 


0.3256 


SSB 

0.2318 

0.2389 

0.2370 

0.2249 

0.2265 

0.2381 

0.2492 

0.2579 

0.2646 

0.2698 

0.2739 

0.2797 

0.2832 

1-5 

A 




0.0278 

0.0280 

0.0295 

0.0310 

0.0324 

0.0335 

0.0344 

0.0351 

0.0362 

0.0368 


LB 



0.0306 

0.0301 



0.0325 





0.0378 



AK 




0.0307 

0.0304 

0.0310 

0.0321 

0.0332 

0.0342 

0.0351 

0.0358 

0.0370 

0.0378 


P 

0.0597 

0.0535 

0.0496 


0.0409 


0.0407 




0.0430 


0.0442 


SSB 

0.0299 

0.0298 

0.0296 

0.0290 

0.0295 

0.0312 

0.0327 

0.0339 

0.0349 

0.0357 

0.0363 

0.0371 

0.0375 

2-3 

A 




1.1925 

1.2239 

1.2592 

1.2884 

1.3107 

1.3275 

1.3404 

1.3510 

1.3679 

1.3821 


LB 



1.1730 

1.2030 



1.2910 





1.3730 



AK 




1.1680 

1.1870 

1.2100 

1.2320 

1.2490 

1.2620 

1.2730 

1.2820 

1.2980 

1.3150 


P 

1.0360 

1.0320 

1.0720 


1.2100 


1.3300 




1.4510 


1.4990 


SSB 

1.1121 

1.0928 

1.1218 

1.1557 

1.1873 

1.2221 

1.2526 

1.2763 

1.2943 

1.3082 

1.3194 

1.3374 

1.3518 

2-4 

A 




0.6848 

0.6785 

0.7010 

0.7279 

0.7518 

0.7716 

0.7879 

0.8014 

0.8221 

0.8368 


LB 



0.7411 

0.7139 



0.7631 





0.8617 



AK 




0.6780 

0.6795 

0.7029 

0.7302 

0.7545 

0.7746 

0.7911 

0.8049 

0.8253 

0.8397 


P 

0.5903 

0.6285 

0.6483 


0.7067 


0.8108 




0.9313 


0.9802 


SSB 

0.6975 

0.7187 

0.7132 

0.6772 

0.6823 

0.7175 

0.7506 

0.7768 

0.7969 

0.8125 

0.8247 

0.8421 

0.8527 

2-5 

A 




0.0833 

0.0840 

0.0884 

0.0931 

0.0972 

0.1006 

0.1033 

0.1054 

0.1085 

0.1105 


LB 



0.0918 

0.0904 



0.0975 





0.1135 



AK 




0.0921 

0.0911 

0.0929 

0.0962 

0.0995 

0.1025 

0.1052 

0.1074 

0.1109 

0.1135 


P 

0.1765 

0.1590 

0.1477 


0.1228 


0.1223 




0.1294 


0.1332 


SSB 

0.0900 

0.0897 

0.0892 

0.0873 

0.0890 

0.0939 

0.0985 

0.1022 

0.1052 

0.1075 

0.1093 

0.1118 

0.1131 

3-4 

A 




1.1413 

1.1308 

1.1683 

1.2131 

1.2529 

1.2860 

1.3132 

1.3357 

1.3702 

1.3947 


LB 



1.2352 

1.1898 



1.2718 





1.4362 



AK 




1.1300 

1.1325 

1.1715 

1.2170 

1.2575 

1.2910 

1.3185 

1.3415 

1.3755 

1.3995 


P 

0.9934 

1.0560 

1.0890 


1.1880 


1.3630 




1.5640 


1.6450 


SSB 

1.1702 

1.2057 

1.1965 

1.1374 

1.1474 

1.2066 

1.2620 

1.3055 

1.3389 

1.3647 

1.3850 

1.4137 

1.4310 

3-5 

A 




0.1388 

0.1401 

0.1473 

0.1552 

0.1620 

0.1676 

0.1721 

0.1757 

0.1809 

0.1842 


LB 



0.1530 

0.1507 



0.1625 





0.1892 



AK 




0.1536 

0.1518 

0.1549 

0.1603 

0.1659 

0.1709 

0.1753 

0.1790 

0.1849 

0.1891 


P 

0.2850 

0.2587 

0.2421 


0.2045 


0.2046 




0.2170 


0.2235 


SSB 

0.1512 

0.1506 

0.1497 

0.1467 

0.1496 

0.1579 

0.1657 

0.1720 

0.1769 

0.1808 

0.1839 

0.1881 

0.1902 

4-5 

A 




0.4708 

0.5463 

0.6114 

0.6468 

0.6630 

0.6687 

0.6692 

0.6670 

0.6599 

0.6524 


LB 



0.4241 

0.4357 



0.5815 





0.6090 



AK 




0.3907 

0.4312 

0.4836 

0.5227 

0.5478 

0.5629 

0.5719 

0.5769 

0.5809 

0.5812 


P 

0.3900 

0.3899 

0.3899 


0.4544 


0.5661 




0.6230 


0.6219 


SSB 

0.3827 

0.3856 

0.3897 

0.4196 

0.5208 

0.5882 

0.6174 

0.6266 

0.6265 

0.6223 

0.6163 

0.6026 

0.5886 
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